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Abstract 

We study quarkonium spectral functions at high temperatures using potential model with com- 
plex potential. The real part of the potential is constrained by the lattice QCD data on static 
quark anti-quark correlation functions, while the imaginary part of the potential is taken from 
perturbative calculations. We find that the imaginary part of the potential has significant effect 
on quarkonium spectral functions, in particular, it leas to the dissolution of the IS charmonium 
and excited bottomonium states at temperatures about 250MeV and melting of the ground state 
bottomonium at temperatures slightly above 450MeV. 
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I. INTRODUCTION 



It is known from lattice QCD that strongly interacting matter undergoes a transition to 
a deconfined phase also called quark gluon plasma (QGP) characterized by chiral symmetry 
restoration and color screening [1]. Quarkonia have been suggested as experimental signature 
of deconfinement by Matsui and Satz |2|. Namely, it has been argued that color screening in 
a deconfined QCD medium will suppress the existence of quarkonium states, signaling the 
formation of QGP in heavy-ion collisions. Although this idea was proposed a long time ago, 
attempts to calculate quarkonium spectral functions from first principle QCD have been 
made only relatively recently. The idea behind these attempts is to calculate the correlation 
function of the corresponding meson operators in Euclidean time in lattice QCD and to 
reconstruct the quarkonium spectral functions using the Maximum Entropy Method (MEM) 
[si- 9] . It turned out that the spectral functions of S-wave quarkonium in the deconfined phase 
are similar to those in the confined phase, and this fact was interpreted as survival of the 
ground state quarkonium in the high temperature region. However, as this was pointed out 
in Ref. ^ details of the spectral functions cannot be resolved at high temperatures due to 
limited extent of the Euclidean time direction (see Ref. [lo[ for recent review on this issue). 
Furthermore, lattice artifacts at small separations (large momenta) further complicate the 



extraction of the spectral functions On the other hand, the temperature dependence of 
the quarkonium correlation functions can be studied reliably in lattice QCD. It was found 
that in the pseudo-scalar and vector channels corresponding to ground state quarkonia 
the temperature dependence of quarkonium correlation functions is either very small or 
moderate [sl, [s], [of. At the same time large temperature dependence in the scalar and axial- 
vector channels was observed. These findings were interpreted as survival of the ground state 
quarkonium and the dissolution of the excited P-states, which fitted well into the expected 
sequential melting pattern. However, a detailed study of the temperature dependence of 
the quarkonium correlators revealed that almost the entire temperature dependence of the 
quarkonium correlators is due to the zero modes not to the dissolution of the bound states 
121. thus challenging the above interpretation. 



With progress in lattice calculations of the correlation functions of static quark anti- 
quark pairs and better understanding of the color screening phenomenon there has been a 
renewed interest in potential models at finite temperature [l3-19|. The basic idea of this 
approach is to assume that medium effects on quarkonium properties can be incorporated 
in a temperature dependent potential and use lattice QCD to constrain its form. It is not 
clear whether and under which circumstances medium effects can be characterized by a 
temperature dependent potential. Fortunately, the effective field theory approach to heavy 
quark bound states provides an answer j20|. Heavy quark bound states at zero temperature 
are characterized by three distinct energy scales m ^ mv ^ mv"^ corresponding to the heavy 
quark mass, the typical momenta inside the bound state, and the typical binding energy 
respectively. The heavy quark velocity v can be treated as a small expansion parameter 
and furthermore, in the weak coupling regime v ~ a^. This allows construction of sequence 
of the effective field theories called NRQCD and pNRQCD by integrating out the scales m 
and mv [2l[ (see Ref. 22 1 for a detailed review on this subject). The degrees of freedom in 
pNRQCD are the singlet and the octet meson fields composed of the heavy quark and anti- 
quark and interacting with ultra-soft gluons, i.e. gluons on energy scale mv"^. The singlet 
and octet potentials appear as the parameters of the effective field theory Lagrangian and 
therefore can be defined at any order in perturbation theory or even non-perturbatively 
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2l| . The effective field theory approach can be generahzed to non-zero temperature, but 
the presence of additional thermal scales T, ~ gT and g^T ( with g being the gauge 
coupling, g^ = Airag) make the analysis more complicated [20[. In the case when the binding 
energy is smaller than the temperature scales the potential receives thermal corrections which 
have both real and imaginary parts and one gets different versions of thermal pNRQCD. 
The precise form of the real and imaginary parts of the thermal corrections depends on the 
relation of the scales 1/r ~ mv, T and gT [20|.^ 



In the past quarkonium spectral functions in QGP have been calculated in Refs. [15 
using potential model with real potentials. The aim of the present paper is to 
perform a calculation of quarkonium spectral functions in the deconfined phase using a 
potential model inspired by thermal pNRQCD, i.e. using a complex potential. The rest 
of this paper is organized as follows: In section [TTl we review what is known about the 
potential at finite temperature from lattice QCD and perturbation theory. In section IIIII we 
present our numerical results for the quarkonium spectral functions. In section |IV] we discuss 
the temperature dependence of the corresponding Euclidean time quarkonium correlation 
functions. Finally, section |V] contains our conclusions. 



II. THE POTENTIAL AT FINITE TEMPERATURE 




FIG. 1: The real (left) and the imaginary part (right) of the potential used in our analysis. 



We are interested in studying quarkonium spectral functions in the temperature region 
where the medium significantly changes the properties of the bound states, i.e. the bound 
states are close to being dissolved. This means that the binding energies are significantly 
reduced from their vacuum values and eventually vanish at some temperature called the 
point of zero binding. Thus, in this situation the binding energy is the smallest scale in the 
problem and all the medium effects can be incorporated in the potential. If we start from the 
pNRQCD and neglect the octet-singlet interactions the dynamics of the singlet quark anti- 
quark field is given by the free field equation which is the form of the Schrodinger equation 



The presence of the imaginary part in the finite temperature potential was first pointed out in Ref. [23j by 
analyzing Wilson loops in the hard thermal loop approximation. In the effective field theory framework 
it can be shown that this analysis corresponds to calculation of the potential in the distance regime 
1/r ^ gT. 
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(e.g. see discussions in Ref. |22|). Then the calculation of the correlation function of singlet 
fields in real time or equivalently the quarkonium spectral function in non-relativistic limit 
is reduced to calculating the Green function of the Schrodinger equation [24 



+ V(r,T)- E 

m 

a{u,T) = -ImG""(f,r"^,E 

TT 



G'''^{f,r',E,T) = 6%r-r') 
E = uj — 2m. 



lr=r'=0 



(2.1) 



The potential entering the above equation is complex [20|, l23[. Its form is quite complicated 
even in the weakly coupled regime as the form of the thermal corrections depends on the 
relation of the scales 1/r, T, mn ~ gT. Furthermore, in the interesting temperature regime 
the separation of the above scales does not hold and the potential is effected by the non- 
perturbative scales g^T, Aqcd- Therefore we have to rely on the lattice calculations to 
constrain the form of the potential. 

The interaction of a static quark and anti-quark (QQ) at finite temperature is usually 
studied on the lattice in terms of the so-called singlet correlation functions Gi{r,T) = 
Tt (W (r)WHO)) with W being the temporal Wilson line. It has been calculated in pure 
gauge theory [isl, 26|, in 3-fiavor QCD [28|, in 2-flavor QCD [2^ and more recently also 
in 2-1-1 flavor QCDfio'-'si'] . Since the expression for Gi{r,T) is not gauge invariant the 
calculations are performed in Coulomb gauge. The singlet correlator can be viewed as 
the correlation function of color singlet static meson field at Euclidean time separation 
T = 1/T. Therefore it is related to the static energy of a QQ pair at finite temperature (see 



discussion in Ref. pi|), though this relation can be quite complicated. The other type of 
static quark anti-quark correlation function that is calculated on the lattice is the Polyakov 
loop correlator G{r,T) = (Triy(r)Triy''(0)) which is gauge invariant. Unfortunately, its 
temperature dependence is much more complicated, because it is sensitive to the excited 
states [29I ^. We also note the properties of the singlet correlator that are discussed below are 
not specific to the Coulomb gauge. In fact, the singlet correlation function defined through 
periodic Wilson loops at finite temperature shares all the properties of the Coulomb gauge 
singlet correlator 3y . 

The logarithm of the singlet correlator defines the so-called singlet free energy 
Fi{r, T)/T = — lnG'i(r, T). For the further discussion it is useful to briefly discuss the prop- 
erties of the singlet free energy. At short distances the singlet free energy coincides with the 
zero temperature potential. In the deconfined region, medium effects become significant at 
distances r,, = O.Aim/{T/Tdec) with T^ec being the deconfinement temperature [25|. At dis- 
tances Te = rT > 0.8 the singlet free energy is exponentially screened, i.e. Fi{r, T) — Foo(T) 
falls off exponentially [2^. The constant Foo(T) is twice the free energy of an isolated static 
quark and it is related to the gauge invariant renormalized Polyakov loop expectation value 
LreniT) = exp {— F^o (T) / {2T)) . We expect that the real part of the potential Rel^(r, T) 



In leading order of perturbation theory the Polyakov loop correlator can be written as the thermal av- 



erage of color singlet and octet contribution [32|, |33| . Recently this has been proved at next-to-leading 



order 



Ref. 



341 ■ where it was also shown how to define rigorously the singlet and octet contribution 



using pNRQCD. In the low temperature limit it can be shown that all excited states contribute to the 
Polyakov loop correlator with the same weight as the ground state 35|, while for the singlet correlator 
the contribution of the excited states is suppressed ^] . 
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entering in Eq. (12. ip will share the same qualitative features as the singlet free energy, in 
particular medium effects will be significant at some distance < rmed{T) < r^. Its value 
will depend on the choice of rmed(T). In our analysis we choose TmediT) = Ve and assume 
that for distances smaller than r^ed the real part of the potential coincides with the zero 
temperature potential, while for r > rmed{T) it is exponentially screened ^. Furthermore, 



we use a smooth interpolation between these two regimes as described in Ref. [19[. This 
procedure also determines the asymptotic value of the potential Voo(T) which turns out to 
be close to value of the internal energy of static QQ pair at infinite separation Uoo{T) (l9| . 
In Fig. [T]we show the real part of the potential for several values of the temperature Since 
on general grounds we expect Fqo < Voo(T) < Uoo the above choice of the potential provides 
the most binding potential compatible with the lattice data. Following Ref. [lH we call this 
choice of Rel^(r, T) the maximally binding potential (in Ref. [19|] it was also labeled as set 
II). Unfortunately lattice QCD cannot yet determine the imaginary part of the potential. 
Therefore here we have to rely on perturbation theory. However, even in perturbation theory 
the form of the imaginary part is only know for certain limiting cases corresponding to some 
hierarchy of the relevant scales 1/r, T, m/j [20[. Therefore for the irnaginary part we choose 
the perturbative hard thermal loop (HTL) result obtained in Ref. |23| 



W^^^(r,T) = - ^ ; " (/>(m^r), 0(x) 



^ stn{zx) 



zx 



(2.2) 



which is in principle valid only for distances r > l/rrio [20|. However, the above expression 
has the nice feature that it vanishes for distances which are much smaller than the inverse 
temperature. For the numerical values of g'^{T) and mD{T) entering the above formulas 



following Ref. [2^ we choose 



Due to the large numerical pref-actor in the argument of the logarithm the imaginary part of 
the potential is not too large even at temperatures T ~ 200MeV. Therefore we may consider 
this choice as a conservative lower bound for it. The resulting Im V"(r,T) is shown in Fig. 
[H In the next section we present charmonium and bottomonium spectral functions for this 
choice of the potential which is maximally binding and minimally dissipative. 



III. NUMERICAL RESULTS FOR THE SPECTRAL FUNCTIONS 

With the potential discussed in the previous section we have solved the Schrodinger 
equation for the non-relativistic Green functions and thus calculated the spectral function 
for S-wave quarkonium. We used the numerical algorithm described in Ref. ^]. First for a 
reference we need to calculate the zero temperature spectral function. At zero temperature 
the potential has no imaginary part, however, Eq. (1) cannot be solved numerically if the 
imaginary part of the potential is strictly zero. Therefore we need to introduce a small 
imaginary part to the potential when solving the Schrodinger equation. In the past we used 
a constant imaginary part equal to 0.01m for charmonium and 0.003m for bottomonium 



Note that the nature of this screening is strongly non-perturbative |37l . [3 
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FIG. 2: The charmonium spectral functions function calculated in potential model neglecting the 
imaginary part of the potential (left) and using both the real and imaginary part of the potential 
(right). The maximally binding potential has been used for the real part. 
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Here we use lmV^^^{r,T = 148MeV)/10 for the imaginary part. It is better 
than the previous choice as it gives less artificial width for the ground state. Furthermore, 
we would like to incorporate the physics of the string breaking into the zero temperature 
spectral functions. Following Ref. [l9| we do this by replacing the Coulomb plus linear form 
of the potential by a Yukawa form at r > l.lfm. The parameters of the zero temperature 
potential can be found in [l9|. For charm and bottom quark mass we choose = 1.29GeV 
and rrif, = 4.67GeV respectively. This gives a reasonably good description of the quarkonium 
spectrum at zero temperature as shown in Table [H The T = spectral functions obtained 
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TABLE I: Spin-averaged charmonium and bottomonium masses in GeV at zero temperature calcu- 
lated in our model and compared to the experimental values f39^. To obtain the 2S spin averaged 
bottomonium mass we used the lattice QCD value of r/;,(25) mass from Ref. {iT 



using the above procedure are shown in Fig. [2] and Fig. [3] for charmonium and bottomonium 
respectively. 

We would like to understand the effect of color screening, encoded in the real part of the 
potential, as well as the effects of dissipation, encoded in the imaginary part of the potential 
on the quarkonium spectral functions. Therefore we first calculated the quarkonium spectral 
functions neglecting the imaginary part of the potential. In terms of the numerics this means 
that we used lmV^^^{r,T)/10 in our calculations. The corresponding charmonium and 
bottomonium spectral functions are shown in the left panels of Fig. [2]and Fig. [3] respectively. 
We see clear peak like structures in the spectral functions which persist to temperature of 
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326MeV for charmonium and 449MeV for bottomonium. In the case of bottomonium we 
also see a remnant of the 25" state in the spectral function at T = 245MeV. This is consistent 



with the results obtained in potential models, in particular, with the calculation of Ref. [17 
that used internal energy as a potential. 

The situation changes dramatically when the imaginary part of the potential ImV^(r, T) = 
lmV^^^{r, T) is taken into account. The corresponding spectral functions are also shown in 
the right panels of Fig. HJand Fig. 12] for charmonium and bottomonium respectively. In the 
charmonium spectral function we see a very broad peak at T = 245MeV, which disappears at 
higher temperatures. In the bottomonium spectral function the 2S peak disappears slightly 
above T = 245MeV, while the ground state peak becomes very broad at T = 450MeV. These 
findings for the spectral functions are consistent with upper bounds on the quarkonium 
dissociation temperatures obtained in Ref. [191] by comparing the temperature dependence 
of the binding energies and the thermal widths, since the lattice data used in both studies 
correspond to a deconfinement temperature of T^ec — ^ 204MeV ^. 

We note that the choice of the real part of the potential used here gives the largest possible 
binding energies. Nevertheless, the binding energies of all quarkonium states, except for IS 
bottomonium are smaller than the temperature for T > 245MeV and drop rapidly with 



increasing temperature [19[. Thus, our approach of putting all medium effects into the 
potential should be at least marginally correct. Binding energies are significantly smaller for 
other choices of the potential. In particular, a theoretically rnotivated choice of the potential 
which is close to the free energy was considered in Ref. [19|, where it was called set I. For 
this potential we see melting of charmonium ground state at T ~ 250 MeV even when no 
imaginary part is considered. At the same time the bottomonium spectral functions are not 
very different compared to the ones calculated with the maximally binding potential. It is 
interesting to see how much our results change when switching the potential to set I and 
keeping the imaginary part of the potential. The corresponding numerical data are shown in 
Fig. m As one can see from the figure changing the real part of the potential has little effect 
on the bottomonium spectral functions when the imaginary part is present. The dissolution 
of the excited states and the broadening of the ground state peak are mostly determined by 
the imaginary part. This also means that the microscopic mechanism behind the melting of 
charmonium and bottomonium states could be quite different. For ground state charmonium 
screening effect play and important role, while for the ground state bottomonium the melting 
will occur due to the large imaginary part of the potential at high temperatures. 



IV. EUCLIDEAN TIME CORRELATION FUNCTIONS 

From the calculated quarkonium spectral functions it is easy to calculate the correspond- 
ing Euclidean time correlation functions 

G(r.T) = r ,.ai.,Tmr.-.n K^r.^.T) = ^"^"'"'^ -//'f' >l (4.1) 

Jo smh(a;/(2T)) 

It is important to compare the temperature dependence of the correlation functions ob- 
tained in potential model with the results of the lattice QCD calculations. The tempera- 
ture dependence of the correlation function on the lattice is studied in terms of the ratio 



^ In Ref. ^19.] the dissociation temperatures were given in units of T^,. 
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FIG. 3: The bottomonium spectral functions function calculated in potential model neglecting the 
imaginary part of the potential (left) and using both the real and imaginary part of the potential 
(right). The maximally binding potential has been used for the real part. 
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FIG. 4: The bottomonium spectral functions function calculated using the complex potential with 
real part corresponding to set I of Ref. 19l |. 



G(T,T)IGrec[j,T), where 



dua{co, T = 0)K{t, u, T). 



(4.2) 



In the lattice calculations this ratio for pseudo-scalar channel that corresponds to S-wave 
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Potential model calculation with real 
Therefore it is interesting to see how close 



quarkonia, was found to be close to one 

potential could reproduce this feature, inereiore it is interesting to see now close is 
G{T,T)/Greci'^,T) to Unity when the imaginary part of the potential is taken into account. 
The spectral functions calculated in the potential model are not reliable far away from the 



threshold because of relativistic effects. Therefore in Refs. |18l . Il9| the spectral functions 
calculated in the potential model were smoothly matched to the perturbative result for the 
spectral functions at some higher energy and the ratios G(r, T)/Grecij, T) have been calcu- 
lated from these matched spectral functions. In this study we did not follow this procedure 
and only used the spectral functions calculated in the potential model when constructing 
G/ (r, T) / GreciT, T). Our results for the ratio of the correlators is shown in Fig. [51 In the case 
of the charmonium this ratio stays close to unity even when the complex potential is used. 
However, G{t,T) /Grec{'T,T) is noticeable smaller compared to the calculations where the 
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FIG. 5: The ratio G{T,T)/GreciT,T) of charmonium (left) and bottomonium (right) correlators for 
different temperatures calculated with complex potential. Open symbols in the left plot correspond 
to the calculation where the imaginary part of the potential was neglected. 



imaginary part is neglected shown as open symbols in the figure. In the bottomonium case 
the imaginary part of the potential had a significant effect and the ratio G{t, T)/Grec{T, T) is 
no longer close to unity in odds with the available lattice data 0, @]- It is possible that this is 
due to fact that time scale relevant for the ground state bottomonium is smaller than the time 
scales related to the medium and modeling of the medium effects by a T-dependent potential 
is not accurate enough. However, the in-medium potential based description should work 
better for the excited bottomonium states due to their smaller binding energies. Therefore 
we constructed the ratio G{T,T)/Grec{T,T) assuming that the ground state bottomonoium 
is unmodified and using only the spectral function above u = 9.7GeV in the calculation of 
the correlation functions. The corresponding results are also shown in Fig. O As one can see 
from the figure the in- medium changes of the bottomonium spectral functions above 9.7GeV, 
in particular, the meting of excited bottomonium states do not produce any large change 
in the Euclidean correlation functions. The problem of bottomonium spectral functions, 
however, needs further investigations in light of recent attempts to calculate t he p roperties 



of ground state bottomonium at high temperatures using perturbative QCD 41 1 and new 
lattice calculations bottomonium correlators in NRQCD at finite temperature 42|. We also 
note that since we do not match the spectral functions to the perturbative form at higher 
energies the calculated G(r, T)/Gj.ec(T, T) shows larger deviation from unity, especially at 



shorter distance compared to previous studies [18|, [1 



V. CONCLUSIONS 



In summary, we have calculated the quarkonium spectral functions in the deconfined 
region using pNRQCD inspired potential model with complex potential. For the real part 
of the potential we considered the lattice QCD motivated maximum binding potential first 
introduced in Ref. [l^, while for the imaginary part we took the perturbative estimate 
from [23|. We find in particular that J/ip peak disappear for temperatures above 245MeV 
and T(IS') peak disappears for temperatures above 450MeV. We found that dissipative 
effects encoded in the imaginary part of the potential play a crucial role in the quarkonioum 
spectral functions and dissolution of quarkonium states. Similar conclusion have been made 
very recently in the analysis of charmonium spectral functions using the T-matrix approach, 
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where dissipative effects are encoded in the heavy quark self energy This study also finds 
no bound state peaks in QGP. In the future it would be interesting to study other choices for 
the imaginary part of the potential, which could lower the dissociation temperature for the 
ground state bottomonium as well as to study the corresponding Euclidean time correlation 
functions more in detail. 
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